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Abstract 

The distribution function of electrons and phonons interacting with electrons in 
semiconductors and semimetals in high electric and quantizing magnetic fields as a 
result of the solution of the coupled system of equations for the density matrixes of 
electrons and phonons is obtained. The effects of heating of electrons and phonons and 
their arbitrary mutual drag are taken into account. The dispersion relation of energy 
of electrons is assumed to be arbitrarily spherically symmetric one. The spectrum of 
phonons is assumed to be isotropic. 

The distribution function of electrons, phonons, the amplification coefficient of 
phonons and the dependence of chemical potential on E, H, n and T e are obtained. 

The distribution function of phonons is obtained for arbitrary drift velocities of 
phonons. 



* e-mail: gassymt@newton.physics.metu.edu.tr 



1 Introduction 



In the present paper, the behavior of semiconductors and semimetals in crossed high ex- 
ternal electric E (E = E x ) and magnetic H (H = H z ) fields is considered. The magnetic 
field is assumed to be high such that the cyclotron frequency Q is much greater than the 
momentum relaxation frequency of electrons v, i.e. 1] C i/. It is known that an elec- 
tron has stationary states at high electric and magnetic fields. In Landau representation, 
the stationary states of electrons are characterized by the magnetic quantum number N, 
the projection of momentum on the magnetic field direction p z and the rotation center of 
electrons X. 

If the spectrum of electrons is assumed to be isotropic and quadratic, then the eigen- 
values of energy have the form[|lJ: 

e a = s N>Pz>x = hn (N + 1) + |S2- - eEX a + . (1) 

The stationary state of electrons with energy given by Eq. (1) is characterized with 

cE 

Hall drift along the y direction with velocity v yo = Vjj = —. In the absence of scattering, 

the average velocity of electrons along the x direction is equal to zero. The inclusion of 

scattering leads to the appearance of the conductance current. The motion of electrons 

along the x direction changes the equilibrium position of the center of oscillation Xq, and, 

as a result, changes the p y component of the electron momentum connected with them 
op 

Xo = --^L. In the case of scattering of electrons by phonons, the motion of electrons 
en 

along the x direction is connected with the transfer of the y component of the momentum 
of electron p y to phonons, and leads to the stream of phonons along the x axis, i.e., it 
leads to the mutual drag of electrons and phonons. During the motion along the x axis, 

cER^n 

the electron gains energy from the field, eE(X a - X,) = eEX a , = — ^ and makes a 

transition from state a to state j3. If this energy is more than the emitted phonon energy 
huj q , electrons and phonons are heated. Thus, the presence of external electric field leads to 
the mutual drag of electrons and phonons. If the electric field is high, then we have heating 
of electrons and phonons and their mutual drag. 

In the quantized magnetic field the scattering frequencies of both electrons by phonons 
Up and phonons by electrons (3 e are increased, (u p , (3 e ~ H 2 ^j . On the other hand, the inter- 
electronic collision frequency v ee is decreased sharply. Because, for nondegenerate statistics 
of electrons only the lowest level of energy is fully occupied, and the number of electrons 
in higher energy levels is exponentially small. In this case the collision frequency between 
electrons in different Landau levels becomes exponentially small. The collisions between 
electrons in the lowest Landau level becomes elastic as a result of the one dimensional- 
ity of the motion, and do not contribute to the interelectronic relaxation. Then, we have 

v eP /T e \ 2 / Ml\ 

— — exp —-^r ■ The decrease of interelectronic collisions in quantized magnetic 



'ee 



Up \hVtJ \ T e J 
fields leads to a sharp decrease of the efficiency of redistribution of electrons in quantum 
states. Therefore, at high E 1 H fields under the reasonable concentrations of electrons, 
the approximation of "effective electron temperature" is not satisfied. That is a reason 



1 



of why for the definition of distribution functions of electrons and phonons it is necessary 
to solve the coupled system of equations for the density matrix of electrons and phonons 
directly. 

Such a problem was formulated and solved for the first time in in the equilibrium 
state of phonons at lattice temperature T. This situation can be realized at high lattice 
temperature T, when the phonon-phonon collision frequency (3 p is much higher than the 
collision frequency of phonons by electrons (3 e . However, for the quantization of orbital 
motion it is necessary to carry out the experiments at low temperatures (liquid Helium, or 
lower) of lattice when (3 P — > 0. As it is shown experimentally0]-0, under these conditions 
it is necessary to take into account the heating and mutual drag of electrons and phonons, 
and the generation of phonons by hot electrons. In the experimental conditions the results 
of satisfy for drift velocities V -C s. Nevertheless, as it follows from the result of 0, 
under the condition V <C s, the heating of electrons is negligible. Moreover, in [Q] the 
change of chemical potential of electrons at high electric and magnetic fields is neglected. 
In fact, as it is shown in the present paper, for the nondegenerate statistics of electrons, 
the presence of high electric field leads to the renormalization of chemical potential. 

The present paper is devoted to solving the coupled system of equations for the diagonal 
parts of electrons and phonons density matrixes at high external crossed E and H fields 
with taking into account the heating of electrons and phonons, and their mutual drag. The 
problem is solved for arbitrary spherically symmetric spectrum of electrons. The spectrum 
of phonons is assumed to be isotropic. 

2 The spectrum of electrons in high electric and quan- 
tizing magnetic fields 

Let us assume that in the absence of electric field, the dispersion relation of electrons is 

B(e a0 ) = ^- + e N , e N = m(N+±y (2) 

The energy of electrons e a at E _L H fields may also be written as 

m(e a )v 2 ya m(e a )m n c 2 E 2 

£ a — £ao — eJ^A a H — e a0 — ei^A a H — — , {6) 

2 m n 2 H z 



where m„, is the effective mass of electrons at the bottom of the conduction band, and 

B(e 

de a 



I dB(e )\ 

m(e n ) = m„ — - — — is the effective mass of electron. In the absence of the electric field, 



Xa = (4) 

m n \l 

For the determination of the dispersion relation of electrons in E _L H fields, we make 

a transition to the reference frame which drifts together with electrons with a velocity of 
cE 

V = — . It is obvious that in such reference frame, the dispersion relation of electrons 
H 

must have the form of Eq. (2), and all properties of the system must be preserved if we 
substitute 
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Pya =P ya + m(e a )V. 
Then, by using Eq. (5) in Eq. (4), we may obtain X a as 



p ya m(e a )V _ p ya eE ( dB(e c 



rn 



,9? \ de r 



(5) 



(6) 



If we define e a o as in Eq. (3) and substitute it into Eq. (2), we may obtain the dispersion 
relation of electrons in E _L H fields as 



B(e*) = B[e a + eEX a 



m(e a ) m n c 2 E 2 



m. 



or, 



2m n 



2 H 2 J 
B(e* a )-e N . 



Pza 

2m„ 



+ hn[N + - 



X a in Eq. (7) is determined from Eq. (6). 

If the increasing of energy of electrons in electric field eEX a — 



m(e n )c 2 E 2 



2 E 2 / 

less than e a , by expanding B(e* a ) into a series around this small parameter, we find 



(7) 
(8) 
is much 



PL 
2m n 



B(e a ) — En + 



eEXa - 



'dBje) 



' dB(e 



or, 



Pla 

2m n 



In the condition V 



m(e a ) m{e a ) m n c 2 E 2 

t>{£a) — £jv H etiX a 7^~TF?- 

m„ m„ 2 H z 



(9) 



(10) 



cE 
IT 



v cr = s, Eq. (10) reduces to the relations obtained earlier 



in [|7| and as a result of the solution of Schrodinger's equation for the Kane spectrum of 
electrons. In the present paper, Eq. (8) is obtained for arbitrary magnitude of the external 
electric field without any limit, and that is why this expression describes a more general 
case than the expressions given in [|7| and [|J. 



3 The main equations and their solutions 

At high magnetic field (Q ^> v) in the Landau representation, the diagonal elements of the 
density matrix of electrons f a is larger than the nondiagonal elements I — ) times, and that 

is why it is enough to write and solve the equation for the diagonal elements of the density 
matrix of electrons. 

Since in the present paper we consider that the space is uniform and the temperature 
and concentration gradients are absent, it is enough to write the equations for the diagonal 
elements of density matrix iV(q, t) for the phonon system. 

After averaging over the electron states with fixed energies, the system of equations 
for the diagonal elements of the density matrixes of electrons f a = f(e a ,t), and phonons 
iV q (t) = iV(q,t) have the form[[|: 
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df(e,t) 2vr 



dt 



X; I C q I 2 I (a I exp(-iq.r) \ (3)\ 2 {5(e p - e a - hu*)[(fp - f a )N(q,t)+ (11) 



fa)} - 6(e p -e a + - U) N M + 



5(£a-£) + J ee [/]+/e d [/], 



^ = yE |C q | 2 |(a|expHqr)|/3)| 2 ^-5 Q -^) (12) 
[(//9 - fa)N(q,t) + fp{\ - f a )\ + I pp [N(q)} + I pb [N(q)}, 

where f(e,t) = X /( e «i t)$( £ a — e) is the average number of carriers with energy e, e.g., 

their distribution function, huf = Tuo q — eE(Xp — X a ), I ee [f] and I e d[f) are the inter- 
electronic and electron-defects collision integrals, respectively, 

IppWJ = P(q)[N(q,t) - N(q,T p )}, J p6 [iV q ] = (3 b [N(q,t) - N(q,T)}, (13) 

are the phonon-phonon and phonon-crystal boundries collision integrals, T p is the tem- 
perature of the heated phonons, and T is the lattice temperature. Both phonon collision 
integrals are considered in the relaxation time approximation: = T p (q), (3 b ~ l = r&fllOl. 



We consider the low electron concentration case n < n cr when I ee <C I ep , where I ep is 
the collisions integral of electrons with phonons. For simplicity we neglect the contribution 
of collisions of electrons with defects. Since in the Born approximation under quantizing 
magnetic field the collision frequency of electrons with neutral and ionized defects do not 
depend on the electron energy, the role of the scattering of electrons by defects may be 
easily taken into account in the final expressions. 

Replacing a by (3 in the second component of the expression under the sum in Eq. (11), 
we find 



dJ ^T L = ¥ £ I I 2 \(a | expHq.r) | /3)| 2 (14) 

{6(ep -e a - K) [Up ~ UMq, t) + fp(l- fa)}} {5(e a - e) - 6(e p - e)} . 

We consider the case when the scattering of electrons by phonons is quasielastic and, 
therefore, changing of the energy of electrons e a — e@ = Tiuj q — Vq y = fko* is less than the 
energy scale of the changing electrons distribution function. 

Then, expanding (fp — f a ) and [6(e a — e) — 8(ep — e)} into series, we find 

//»-/« = f{e a + K) - /(*«) « H ( ^4^-) + (15) 



de. 



a 



d 

{5(e a -e)- 6(6/3 - e)} = S(e a -e)-S(e a -e + fwj* q ) « -huj* — 5(e a - e) + (16) 
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d d 
By using the identity — — 5(e a — e) = —5(e a — e), we obtain 

L/C. f\ LyCs 



where, 



= yE I c q I 2 14/3 1 2 (K) 2 ^ ~ e «- k) ^ ( £ « - £ ) > ( 18 ) 

a,/3,q 

= "T" E I C <i I 2 ! J «/3 I 2 ^9 5 ( £ « ~ £ /3 - Huj g) 6 ( £ « - £ ) ■ 

11 a,/3,q 

The stationary solution of Eq. (17) satisfying boundary condition lim^oo f(e) — > is 

/( £ ) = jcon^exp^^^+l} , (19) 

where T e (e) = is the temperature of electrons with energy e. 

The solution of Eq. (12) is 

N(q, t) = {N(q, 0) + (3^ N(q)} exp ( 7 ,i) - /^J 1 iV(g), (20) 

where 7 „ = /3 ( - — — 1 J is the increment of the generation of phonons, (3 = (3 e + (3 V + is 
the total collision frequencies of phonons by the scatterers and 



& = X E I ^ I'! 4/3 I 2 (4 - /a) * (e/3 " e a - K) ( 21 ) 



« f E I C q | 2 | lap I 2 (e/j - e«) f (e/J - e« - K) ■ 

«,/3,q \ a / 

iV(g, 0) = N(q, T) is the initial distribution function of phonons at t = in the absence 
of external fields, 

«=(l-^)v, |/ Q/3 | 2 =|(a|e^Hq.r)|/3)| 2 . (22) 

As it follows from Eq. (20) in the 7g > case, i.e., when the drift velocity of phonons u 
is larger than the sound velocity s, the distribution function of phonons iV(q, t) increases 
exponentially with time, whereas in the 7g < case the solution, Eq. (20), is stationary, 

7V(q) = lim JV(q, i) = -^N^q), iV(g) = 7e iV(q, T e ) + lp N(q, T p ), (23) 
where 7e = — , 7p = — , 
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£ I I aP | 2 T c (e) [9/(e a )/9e a ] ^ - e a - ftu;*) 

^(g,T e ) = ^ , (24) 

2^ | I a/3 | ftw* [df(e a )/de a ] 5{e p -e a - hu*) 

a/3 

I 4/3 I 2 r e (e) [df{e a )/de a ] 5{ep -e a - hcu*) 

T e = - 2 . (25) 

2^ | I aj 3 I [df(e a )/de a ] ^(ep - e a - ftw,) 

a/3 

where T e is the effective temperature of electrons. 

u d f is t\ 

In the case of ^ q = or — — 1 or hu* = 0, from Eq. (17) we have — — - — = 0, or 

TlUlq C/t 

f(e, t) = const. At this point the distribution function of electrons f(e, t), the effective elec- 
tron temperature T e and the drift velocity of electrons are constant. Therefore, the state is 
nondissipative and the current is constant. Moreover, the distribution function of phonons is 
nonstationary and grows by time linearly, N(q, t) = N(q, T) + (3tN(q, T e ). Indeed, u = V = 

s = const, and dN ^ - = (3N(q,T e ) = const. Namely, P(T e ) = ^hw* (^—^^j = °> 

where P(T e ) is the power transferred by electrons to phonons. The point u = s is the 

acoustical instability threshold (AIT). At this point, the stimulated emission of phonons is 

equal to the stimulated absorption of phonons, and we have only spontaneous emission of 

phonons at high external electric and magnetic fields. At this point collisions of electrons 

with phonons are exactly elastic, i.e., state is nondissipative and dynamically stationary 

because of the power received from electric field emitted as phonons by the process of the 

cE 

stimulated emission. Moreover, at this point the current J = const, and v d = V = — = s, 

H 

sH cE 
or cE = sH, i.e., E cr = — or H c 



L cr 

C 



Substituting Tej = in Eq. (12), we may directly solve this equation under the 

(J Is 

boundary conditions 

iV(q,T e )| t=0 =iV(q,0)^iV(q,T) = |exp^^-l| , (26) 

and we obtain Eq. (23). 

Under the considered conditions hu* <C T e ,T from Eq. (23), we find 

~ q) = 7e^+7pT p = T_ = ~ = ^ + ^ (27) 

where T is the temperature of coupled by the mutual drag system of electrons and phonons. 
Therefore, the stationary solution of Eq. (12) has the form (hu* > 0): 



* (q) = ^rq = ~P N (v) « (28) 

flUJa 
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Let us substitute Eq. (28) in Eq. (18) and take into account the relations: 



'a/3 1 



1/2 



l NN' 



N'.\ 



2 - 1 T 
— \-LNN' 


I 2 *, 


-*)(■ 


_Q± 







(29) 



N'-N 



Qh 



Here = hR" 1 = (m^^) 1 / 2 , q± = ql + 0%, -^iv - ' is the Laguerre's polynomial normalized 
to unity, and | m \ = \ N — N' \. 

Choosing the cylindrical coordinate system with direction along the magnetic field 



dq = q± dq± dq z dip, q y = q± sin ip. 
As a result of the integration, we have 



(30) 



A(e) = T<$>{u/s)D{e), <$>{u/s) = <p - 2^ - 1) - - l) 2 , (31) 

x 2 | Inn> | 2 m 2 (e) 



1/2 



B(e + Tllu*) — e N i B{e — huj*) — en 



1 1/2 ' 



where x = q±/qH, and W(qn) is the constant part of the potential in mutual interaction of 
electrons by phonons. 

Thus, in general, under the arbitrary degree of quantization, the ratio of A(e) to D(e) 
is T$(u/s); and it does not depend on the energy of electrons and the potential of mutual 
interactions with phonons. In other words, in a more general case of interaction of electrons 
by the acoustic and optical phonons the expression: 



T e/ = ^ =?*(«/*), (32) 



does not depend on the energy of electrons. 

By substituting Eq. (32) into Eq. (19), we find 



m J 1+ ^(i^l)Y\ (33) 



In other words, at quantizing magnetic fields the distribution function of electrons is a 
Fermi one, and includes the effective temperature in the following form: 

T e/ = f{l+(l-^)V-l)}, ¥>i=(l-£) U2 ■ (34) 



In the classical region of strong magnetic fields (Q u) 
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t(M/s)= { 1+ ( 1 _^) fe _ 1 ) + lQ 2 j, V2= _l ln 

As it follows from Eq. (33), the distribution function of electrons is the Fermi one with 
effective temperature: 



s + u 



u 



(35) 



The fact that the distribution function of electrons in both classical and quantum re- 
gions of magnetic fields are the Fermi ones with effective temperature is a result of the 

cE 

independency of drift velocity of electrons V = — from energy of electrons e. 

H 

At small values of the drift velocity of electrons, V <C s, from Eq. (34) for the quantizing 
magnetic field, we have 



If p e » fa, T ef = f jl - i (^|) 2 }; and if fa « fa, T ef = f jl + i (|§) 2 }. 
In the classic region of magnetic fields 

If fa 3> fa, phonons are not heated, T = T and F C s. Under these conditions from 
Eqs. (37) and (38) we can obtain the results of 0. 



4 The statistics of electrons at high magnetic fields 

In crossed E _L H field, the chemical potential of electrons may be obtained from the 
normalization condition of the distribution function as 

n e = fds g H (e)f(e), g H {e) = ZX^)- (39) 

For arbitrary spherical symmetric spectrum of electrons 



2(2m n ) 1 / 2 m n n _ [°° , 1/2 f 

= — 2^/ [S(e ) - ejv] 7 < 1 + exp 

n Je * I 



(2tt^) 2 



'e*-C(E,HY 



T, 



<-f 



C(E,H) = C(E,H) + eEX(e)~ 



m{e)c 2 (E^ 
2 [lP 



(40) 
(41) 
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The £* is determined as a solution 

B{e*) -e N = 0. 
By partial integration of Eq. (40), we find 

4(2m n ) 1 / 2 m„fi_ /■«> 



_ /"OO 



1/2 



(2tt^) 2 

For the case of parabolic spectrum of electrons 

e\ = e N = TiVL (n + -) , 
and for the Kane spectrum of electrons: 



9/(e*) ' 



1-1 



4e N 



1/2 



(42) 



(43) 



(44) 



(45) 



Let us consider the case of degenerate and non-degenerate statistics of electrons alone. 
For strong degenerate electrons 



' dfoje* 
de* 



S(e - C 



With the help of this expression, we can integrate the Eq. (43) and get 



4(2m n )V 2 
n e = _ tn 2^ 



h(2nR) N 
For the parabolic spectrum of electrons 



1/2 



4(2m ra )V 2 

AT 



ft(27ri?) 2 

For the Kane spectrum of electrons 



1/2 



4(2m„)V 2 

//, - 2^ 

N L 



h{2-KR) 2 



C (l + —) -mi(tv + J 



,1/2 



For ultraquantum limits (N — N f — 0), we have 



(46) 



(47) 



(48) 



(49) 



n P = 



4(2m„ 



|V2 



B(C) 



1 1/2 



_ nn ft 2 (27riZ) 4 n 2 
l ^ J 2 16(2m n ) ' 



(50) 



/i 2 f27ri?) 4 7?, 2 

where B(C*) = 1 — -. From this relation for the parabolic spectrum of elec- 

2 16(2mn) 

trons we obtain the chemical potential as follows: 
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(51) 



This expression is the same of Eqs. (32) and (22) in [12 
In the case of Kane spectrum of electrons: 



C(E,H) = -■£{! 



2hQ 



27r A h 2 R A n 2 \ 1/2 ' 



nine. 



(52) 



In the case of nondegenerate electrons by taking [B(e*) — sn] out of the integral when 
e* = e\(N) + T e , we find 



exp 



'(**(E,H)\ (2nR) 2 hn e 



4(2m n )V2 



exp 



N 



4(N) ' 



[B{e*{N)+T e )-e N 



1/2 



(53) 



As before e\ is obtained from Eqs. (44) and (45). If the condition e\ ^> T e is satisfied, 
then by expanding B(e\ + T e ) in Eq. (53) into series and by taking into account Eq. (42), 
we have 



exp 



'(**{E,H)\ (2nR) 2 hn e 



T 



AV2 



rp-l/2 



exp 



L N 



-i -l 



where 



m{e\) = m n 



'dB{e* 



(54) 
(55) 



Dividing Eq. (53) by itself for E = 0, we get 

T. 



where, 



C(E, H) = C(H)^ + T e In + T e In 

m m n c 2 E 2 



F N (T) 
F N (T e ) 



r(E,H) = ((E,H)- 



m n 2 H 2 ' 



MTe) = £ exp (-^p) {B(el(N) + T e ) - e N } 



1/2 



(56) 

(57) 
(58) 



Fn(T) may be obtained from Eq. (58) by replacing T e with T. 



For the case of parabolic spectrum of electrons e\ = = ^Vt ( iV + - ) , and, we have 



(**{E,H) = C{H) J ± + T e ln 



n 



T 

-^ln 

2 



+ T e ln 



sinh(^fi/2T e 



sinh(M2/2T) 



(59) 



Therefore, in the case of parabolic spectrum of electrons we finally find the chemical po- 
tential as 
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«E, H) = C(H)^ -eEX + '^-^+ T e In 



2 VT 



+ T e ln 



sinh(ftft/2T e ) 



+ 



(60) 



sinh(M2/2T) 



As it follows from Eq. (60) in external electric field if the concentration of electrons is 
increased (n e > n), then the chemical potential of electrons must also increase. 

If we have full ionization of small impurity centers, then n e = n = const. Thus, from 
Eq. (60), we may obtain 



((E,H)- — = ^ 



In 



+ T e ln 



m„c 2 E 2 
_l 5 

2 H 2 
exp(hn/2T e ) 
exp(htt/2T) 



(61) 



In the quantizing magnetic fields hQ > T e ,T and, for this reason, the last term in Eq. 
(61) may be presented as m this case 



~2~ 



(62) 



For the weak electric fields T P 



(63) 



As it follows from Eq. (63), the expression for ((E,H) differs from the expression 

TYl C 2 E 

given in |13| by the factor — - — —, which is connected with the Hall drift of electrons. 



2 H 2 ' 

In the case of heating of electrons at external electric field the main contribution to the 

T e 

free energy is obtained by the expression ((H)—, i.e., in high external electric field the 

T e 

chemical potential of electrons in common case increases linearly by — . 

For the calculation of the statistical behavior of electrons in external fields such as 
magnetic susceptibility, heat capacity and etc., it is necessary to know the dependence of 
the chemical potential on the intensity of external electric and magnetic fields. 

The magnetic susceptibility of hot electrons in the case of high concentration of electrons 
was investigated earlier |TJ[ . 



5 Conclusion 



In the present work, it is shown that under the conditions of arbitrary degree of quanti- 
zation and for the interaction of electrons with both the acoustical and optical phonons, 
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the distribution function of electrons has the form of Fermi distribution function with ef- 
fective electron temperature. This result is obtained in the case of small concentration of 
electrons n < n cr when the usual approximation of "effective temperature" for electrons is 
not satisfied, i.e., u ee <C v v . 

The distribution function of phonons interacting with electrons is obtained for arbitrary 
drift velocities of phonons. If the drift velocity of phonons u is smaller than the sound veloc- 
ity s there is a stationary state and the distribution function of phonons is shifted Planck's 
one with effective temperature of phonons. In the case of u > s, the distribution function of 
phonons grows with time exponentially, i.e., we have the effect of generation or amplification 



of phonons by electric field. The amplification coefficient of phonons 1 q = — \ t 1 is 



obtained. It has been established that in the considered case the chemical potential of 
electrons is renormalised. The E, H, n and T e dependences of the chemical potential are 
obtained. 
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